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The two-fragment electrodisintegration of 4 He into proton 
and triton is calculated in Plane Wave Impulse Approxima- 
tion (PWIA). The three- and four-nucleon wave functions in- 
volved are obtained by solving the Alt-Grassberger-Sandhas 
(AGS) integral equations, with the Malfliet-Tjon potential 
as the underlying NN-interaction. Our results are in re- 
markable agreement with the experimental data and, in con- 
trast to alternative approaches, do not exhibit any dip in 
the five-fold differential cross section at a missing momen- 
tum of ~ 450 MeV/c. 



PACS numbers: 21.45. +v, 25.10.+S, 25.30. Fj 



The two-fragment electrodisintegration process 
4 He(e, e'p) 3 H has been the subject of several experimen- 
tal investigations for various kinematics (see, for exam- 
ple, [|]~^)). On the theoretical side quite some effort has 
been devoted to calculating this process. However, the 
exact treatment of four-nucleon electrodisintegration ob- 
servables is computationally very demanding and, thus, 
has usually been simplified by approximations and model 
assumptions |||-|| . 

In Plane-Wave Impulse Approximation (PWIA) all 
these calculations exhibit a characteristic dip, actually 
zero, in the five-fold differential cross-section around a 
missing momentum of ~ 450 MeV/c, which does not show 
up in the experimental data. Laget performed calcu- 
lations including final state interaction (FSI) effects and 
meson exchange currents (MEC) by means of a Feynman 
diagrammatic approach. Although this resulted in a par- 
tial filling of the dip, these investigations also underesti- 
mate the data considerably in this region. Similar results 
were obtained when FSI was taken into account via an 
effective nucleon-trinucleon interaction [ ppH I . In a com- 
pletely different approach Nagorny et al. |10| ] included the 
electromagnetic field within the strongly interacting sys- 
tem in a relativistic gauge invariant way, the FSI being 
incorporated via the pole contribution of the p 3 H— + p 3 H 
scattering matrix fTl) . The agreement with the data is 
again fairly satisfactory, but the zero is exhibited as well 
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In detail, at missing momenta less than 300MeV/c all 
calculations show a good agreement with the data [§,[|,D . 
Surprisingly, the PWIA performs reasonably well in this 
region where the FSI could be expected to be more im- 



portant than in the higher missing momenta region. In 
contrast, in the region 300MeV/c< Q < 600MeV/c the 
results strongly depend on the way the FSI effects are 
included. For example, as pointed out in || the Laget 
results underestimate the cross section by a factor of 4 
and those of Schiavilla by a factor of 2. At missing mo- 
menta above 600MeV/c, where the MEC contribution 
is becoming important, the agreement with the data is 
again fair. We, therefore, conclude that the zero in the 
PWIA cross section is not necessarily a manifestation of 
strong FSI or MEC effects. 

Instead, one should look for other explanations, such as 
the dependence of the results on the model used, the NN 
forces employed, the determination of the bound state 
wave functions etc. In the present work the wave func- 
tions involved were calculated within the exact three- and 
four-nucleon AGS formalism |I^,[l3|]. We mention that 
the same wave functions have successfully been employed 
already in calculations of the two-fragment photodisinte- 
gration of the a-particle fb|,[L5| . 

We consider the two-fragment reaction in which the 
scattered electron and the ejected nucleon are measured 
in coincidence. The corresponding electron-proton coin- 
cidence cross section is given by 
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where ctm is the Mott differential cross section, 
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Ei(Ef) is the energy of the incoming (outgoing) electron 
and pf is the relativistic density of states. The transi- 
tion matrix, properly antisymmetrized with respect to 
the four nucleons Jig], is given by 



A4(q) = 2<->(q;*/jj|ir|¥/v) 



(3) 



where H is the Hamiltonian describing the interaction 
between the electron and the nucleons. The ejected pro- 
ton moves away with momentum q with respect to the 
residual three-nucleon bound state \^fm)- The kinemat- 
ics of this process is shown in Fig. 1. 
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The Hamiltonian for the interaction between an elec- 
tron and four nucleons is that of McVoy and van Hove 
JT7| , which has been previously employed in the electro- 
disintegration of the trinucleon system by Lehman and 
collaborators and by Epp and Griffy p|. This 

Hamiltonian, correct to the order of h 2 Q 2 /M 2 cr, is 
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Here Xj and are the position and momentum opera- 
tors of the j-th nucleon, <Tj is the nucleon spin operator, 
a is the Dirac matrix acting on the free electron spinors 
\vi) and \vf), while q 2 ^ is the exchanged four- momentum 
squared. Fin and F 2 n are the form factors of the nu- 
cleon, k is the anomalous moment of the nucleon in nu- 
clear magnetons, and M is the nucleon mass. 

For proton knock-out, the transition matrix Eq. (||) 
reads 

M = -(vf\vi)M Q + {v f \a\vi) • (Mei + M mag ) , (5) 
where 

M q = 2^(<i^ iii \Hq\^ iv }, (6) 
M el = 2(-)(q;* /I j|Hei|*jv), (7) 
M mag = 2(-)(q;* m |H mag |* 7V ). (8) 

The Hamiltonians Hq, H e i, and H mag are given by 



n Q = F^(l + ql/8M 
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Here the superscript p refers to the proton and Xj = 
(1 +r|)/2 is the isospin operator for nucleon j while F^ h 
and -F£ ag are the charge and magnetic form factors of 
the proton defined by 



^ch = F iP + (ql/4M 2 ) Kp F 2p 
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The analytical fit to the proton form factors F± p and F 2p 
given by Janssens et al. plf is used in the calculations. 



Squaring the matrix element, summing and averaging 
over the electron spin, and inserting the resulting expres- 
sion in Eq. (pi), we obtain 
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where J = M c i + M mag . The determination of the coinci- 
dence cross section is thus reduced to the determination 
of the nuclear matrix elements A4q and J. 

In this work we use the dominant electric operators (|^) 
and (|h]). In PWIA the nuclear matrix elements (^|) and 
(@) read 

4 

B Q (q) =2<q|(tfjj J |£exp(-iQ-x i )A i |* / v) (15) 
3=1 

and 

4 

Bei(q) = 2 (q|(*///| (Pi exp(-iQ ■ xj) 
3=1 

+ expHQ • Xj)Pi) Aj|*jy) , (16) 

where q = (pi + p 2 + P3 - 3p 4 )/4. 

The operators appearing in Eqs. (|l5|) and (|l^) are 
the same as those of Eqs. (^) and (|l(]), except that the 
nucleonic form factors F^ and i^ag are n0 * noted here. 

To proceed we express the operators Xj and Pj in Ja- 
cobi coordinates and neglect, as in the photodisintegra- 
tion case [14 ll|, those acting within \^fiu). The re- 



maining terms, containing q and its canonically conju- 
gate counterpart, are treated without further approxi- 
mation. A straightforward calculation then reduces jl5| ) 
to 



1 



Bg(q) = 2 <q + -Q|(*///|(Ai + A 2 + A 3 )|* 



4 

2<q-^QK* /ZI |A 4 |*jy) ! 



IV 



(17) 



whereas ([L6|) is replaced by 
4 1 

B e i(q) = gq(q+^Q|(4'///|(Ai+A 2 + A3)|*/y) 

-4q(q-^QK* m |A 4 |*zy>+ QB Q (q). (18) 

The construction of the above matrix elements requires 
the knowledge of the bound states \^m) and l^/y). For 
their calculation the exact three- and four-nucleon AGS 
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integral equations are employed |LJ,|L3j. The latter con- 
sist the coupled set of (before antisymmetrization) 18x18 
four-body AGS equations. They contain in their kernel 
all subsystem information via the two-body T-matrices, 
the three- and (2+2)-body AGS transition operators. By 
this approach, the full coupling and the corresponding 
interference of (2+2)- and (3+l)-channels in the four- 
body system is taken into account explicitly, and thus 
exactly and completely. In order to reduce the origi- 
nal three- and four-body relations to (one-dimensional) 
integral equations, the W-matrix method p3| and the 
energy-dependent pole approximation (EDPA) (23| are 
used. In the purely nuclear case these approximations 
have led to very accurate results (see e.g. |2^-^]). Fur- 
thermore, they have been successfully used in calcula- 
tions of the photodisintegration of 3 H, 3 He [| 27 and 
4 Hc p^ , |l5| ,p| . The graphical representation of matrix 
elements like Eqs. (|l5|) and (|l6|), adapted to the four- 
body AGS formalism, can be found in fill . As in fl4[ , |l5|l , 
the Malfliet-Tjon potential I and III |2q | is chosen, as 
it is both sufficiently realistic and simple enough to be 
employed in four-nucleon calculations. This property is 
of particular importance in our calculations where the 
computation of the matrix elements, despite the approx- 
imations used, is still tedious and of considerable numer- 
ical complexity. The corresponding binding energies are 
8.595 MeV for 3 H and 30.1 MeV for 4 He Q. 

The results obtained within the AGS formalism for 
the 4 He(e,e'p) 3 H five-fold differential cross section as 
a function of the missing momentum Q are shown in 
Fig. 2. The kinematics and the experimental data are 
those of HH for the ui = 215 MeV case. For compari- 
son we also included in the figure the PWIA results of 
obtained for wave functions constructed via the in- 
tegrodifferential equation approach (IDEA) of Ref. [p9| . 
The PWIA results of Laget (see Refs. @,|,@) are also 
shown, being obtained for the Urbana potential and for 
wave functions constructed with the variational Monte 
Carlo (MC) method. The agreement of our AGS cal- 
culations with the experimental data, especially in the 
region where the PWIA results of the other two methods 
show their characteristic dip, is remarkable. This holds 
true also in comparison with the other results reported 
in ||. Fig. 3 shows the five- fold differential cross section 
for the 4 He(e,e'p) 3 H reaction for the Saclay kinematics 
H . For comparison the Laget results are also plotted. 
The agreement of our results with experiment is again 
remarkable. 

The overall small discrepancies may be reduced by us- 
ing a better NN force, further improvements in the PWIA 
matrix elements, and inclusion of the FSI in a rigorous 
way. A particular advantage of the AGS-type approach 
in this respect lies in the fact that the incorporation of 
the FSI is exactly the same for the four-nucleon scatter- 
ing, the photodisintegration of 4 He, and the electrodis- 
integration of 4 He (see e.g. and Refs. therein). Most 



important: in this approach the underlying integral equa- 
tions explicitly incorporate the (2+2)-channels, not fully 
included in other approaches, and their coupling to the 
(3+l)-channels. That underlying interference of compet- 
ing channels is the most relevant feature of four-body the- 
ory as compared to three-body theory. In other words, 
the complexity of four-body rearrangement processes is 
fully taken into account. 

In conclusion, our results show that already in PWIA 
quite a good description of the experimental data can be 
achieved. The main reason for this agreement appears 
to be the use of wave functions obtained from the AGS 
integral equations with their complete coupling scheme. 
Another reason is the way of calculating the nuclear ma- 
trix elements. Namely, those parts of the electromagnetic 
operators @ and ([l0|), which act between the relative 
motion of the two outgoing nuclear fragments, are taken 
into account exactly. The sensitivity to the input NN- 
potential and to the above-mentioned 2+2 rearrangement 
terms are under investigation. 
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FIG. 2. Five-fold 4 He(e,e' p) 3 H differential cross section 
as a function of the missing momentum Q for the NIKHEF 
kinematics [4,5]. 
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FIG. 3. Five- fold 4 He(e,e'p) 3 H differential cross section 
as a function of the missing momentum Q for the the Saclay 
kinematics [3]. 



FIG. 1. Kinematics for the process 4 He(e, e'p) 3 H. 
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